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Context. The magnetic fields of the accreting white dwarfs in magnetic cataclysmic variables (mCVs) determine the accretion geometries, the 
emission properties, and the secular evolution of these objects. 

Aims. We determine the structure of the surface magnetic fields of the white dwarf primaries in magnetic CVs using Zeeman tomography. 
Methods. Our study is based on orbital-phase resolved optical flux and circular polarization spectra of the polars EF Eri, BL Hyi, and CP Tuc 
■~>«^ obtained with FORSl at the ESO VLT. An evolutionary algorithm is used to synthesize best fits to these spectra from an extensive database of 
pre-computed Zeeman spectra. The general approach has been described in previous papers of this series. 

Results. The results achieved with simple geometries as centered or offset dipoles are not satisfactory. Significantly improved fits are obtained 
^ for multipole expansions that are truncated at degree /„,ax = 3 or 5 and include all tesseral and sectoral components with < m < I. The most 
frequent field strengths of 13, 18, and lOMG for EF Eri, BL Hyi, and CP Tuc and the ranges of field strength covered are similar for the dipole 
C/j and multipole models, but only the latter provide access to accreting matter at the right locations on the white dwarf. The results suggest that 
the field geometries of the white dwarfs in short-period mCVs are quite complex with strong contributions from multipoles higher than the 
^ ^ dipole in spite of a typical age of the white dwarfs in CVs in excess of 1 Gyr 
• , Conclusions. It is feasible to derive the surface field structure of an accreting white dwarf from phase-resolved low-state circular spectropo- 
' larimetry of sufficiently high signal-to-noise ratio. The fact that independent information is available on the strength and direction of the field 
^ in the accretion spot from high-state observations helps in unraveling the global field structure. 

Key words. stars:white dwarfs - stars:magnetic fields - stars:atmospheres ~ stars: individual (EF Eri, BL Hyi, CP Tuc) - polarization 



1. Introduction 

The subclass of magnetic cataclysmic variables (mCVs) termed 
polars (Krzeminski & Serkowski 1977) contains an accreting 
white dwarf that emits circularly polarized cyclotron radiation 
from an accretion region standing off the photosphere, often 
refeiTed to as accretion spot. The harmonic structure of the 
cyclotron radiation allows a straightforward measurement of 
the magnetic field and an estimate of the field direction in the 
spot. Photospheric absorption lines are heavily veiled by the 
intense cyclotron emission in the high (accreting) state. The 
field structure over the surface of the white dwarf becomes ac- 
cessible to measurement only in low states of discontinued ac- 
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cretion via the profiles of the photospheric Zeeman-broadened 
Balmer absorption lines, an approach which is applicable also 
to non-accreting isolated white dwarfs. Different from the lat- 
ter, accreting systems offer the advantage that the strength and 
direction of the field in the accretion spot and its approximate 
location on the surface as determined from high-state observa- 
tions represent a fix point for the field structure. For simplicity 
it was often assumed that the field is quasi-dipolar, although 
accreting systems with two accretion spots separated by much 
less than 180° supported suspicions of a more complex struc- 
ture (Meggitt & Wickramasinghe 1989, Piirola et al 1987b, see 
Wickramasinghe & FeiTario 2000 for a review). 

We have set out on a program to obtain a more complete 
picture of the surface field structure of magnetic white dwarfs 
using an approach dubbed Zeeman tomography (Euchner et 
al. 2002). The field geometries of two isolated white dwarfs, 
HE 1045-0908 and PG1015+014 (Euchner et al. 2005, 2006), 
proved to be significantly more complex than simple centered 
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or offset dipoles. In this paper, we present first results of the 
Zeeman tomography of three polars observed in their low 
states, EF Eri, BL Hyi, and CP Tuc, and find that they, too, have 
rather complex field geometries. 



2. Observations and Data Analysis 

We have obtained spin phase-resolved circular spectropo- 
larimetry of EF Eri, BL Hyi, and CP Tuc in their low states. 
These stars belong to the short-period variety with orbital pe- 
riods of Sl.Omin (EFEri), 113.6min (BL Hyi), and 89.0min 
(CP Tuc). The secondary star is a late M-star in BL Hyi. It is 
substellar in EF Eri (Beuermann et al. 2000; Harrison et al. 
2004) and possibly in CP Tuc, too. Full orbital coverage was 
achieved for EF Eri and BL Hyi, but for technical reasons only 
half of the orbital period was covered for CP Tuc. 

The data were collected at the ESO VLT using the focal 
reducer spectrograph FORS 1 . The instrument was operated in 
spectropolarimetric (PMOS) mode, with the GRIS_300Vh-10 
grism and an order separation filter GG 375, yielding a usable 
wavelength range of 3750-8450 A. With a slit width of 1" the 
FWHM spectral resolution was 13 A at 5500 A. A signal-to- 
noise ratio of typically S /N - 100 was reached for the individ- 
ual flux spectra. FORS 1 contains a Wollaston prism for beam 
separation and two superachromatic phase retarder plate mo- 
saics. Since both plates cannot be used simultaneously, only 
the circular polarization has been recorded using the quarter 
wave plate. Spectra of the target star and comparison stars in 
the field have been obtained simultaneously by using the multi- 
object spectroscopy mode of FORSl. This allows us to derive 
individual correction functions for the atmospheric absorption 
losses in the target spectra and to check for remnant instrumen- 
tal polarization. Tab. 1 contains a log of the observations. 

The observational data have been reduced according to 
standard procedures (bias, flat field, night sky subtraction, 
wavelength calibration, atmospheric extinction, flux calibra- 
tion) using the context MOS of the ESO MIDAS package. In 
order to eliminate observational biases caused by Stokes pa- 
rameter crosstalk, the wavelength-dependent degree of circular 
polarization V/I has been computed from two consecutive ex- 
posures recorded with the quarter wave retarder plate rotated 
by +45°. The circular polarization was then obtained as the av- 
erage of two consecutive sets of spectra in the ordinary and the 
extraordinary beams (see Euchner et al. 2005, for details). 

The three stars were in their low states with magnitudes 
estimated from the spectrophotometry ofV~ 18 for EF Eri and 



Table 1. Dates of the spectropolarimetric observations obtained at the 
ESO VLT, exposure times, and number of exposures. 
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Fig. 1. Flux spectra (top) and circular polarization spectra (bottom) of 
EF Eri at two selected orbital phases and for two field models. Model 
A is an off^set dipole and model B a multipole expansion truncated at 
/max = 5. The data are shown as grey curves, the best fit models from 
the Zeeman tomographic analysis are overlaid as solid black curves. 
The order of orbital phases and field models for the polarization spec- 
tra is the same as for the flux spectra, i.e., the bottom spectrum and 
third from the bottom are for the multipole expansion, the other two 
for the ofi'set dipole. The ordinate scales refer to the bottom spectrum, 
the other ones are arbitrarily shifted upwards. 



BL Hyi and y ~ 19 for CP Tuc. CCD photometry of BL Hyi in 
the same night gave V = 1 7 .45 . In case of BL Hyi, the resulting 
spectra were corrected for the contribution by the secondary 
star using a spectrum of the dM5.5 star G1473 as a template. 
No trace of the secondary star was seen in the other two objects. 
Seeing variations and a loss of blue flux in the first two hours of 
the EF Eri run and in the last three spectra of BL Hyi afi'ected 
the detection of orbital modulations. The orbital modulation 
of EF Eri seen in the remainder of the data is consistent with 
that reported by Szkody et al. (2006, and references therein). 
No substantial orbital modulation was detected in the data of 
BL Hyi and CP Tuc. 

In order to facilitate analysis of the Zeeman absorption fea- 
tures, the continua of the observed spectra were normaUzed by 
the following procedure which minimizes the differences be- 
tween observed and theoretical continua and corrects also for 
the mentioned loss of blue Ught. In a first step, a mean ef- 
fective temperature of each object was determined by fitting a 
magnetic model spectrum to the quasi-continuum of the mean 
of the observed spectra unaffected by light loss. In a second 
step, the continua of all observed flux spectra were sampled 
in some 20 narrow fiducial wavelength intervals, which avoid 
the known Zeeman features and the emission lines, and ad- 
justed to the best-fitting model spectrum using low-order poly- 
nomials (see Euchner et al. 2005, for more details). This ap- 
proach removes the time variability in the flux continua but 
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Fig. 2. As Fig. 1, but for BL Hyi and CP Tuc. Model A is an offset dipole and model B a multipole expansion truncated at /max = 3. 



leaves the equivalent widths of the Zeeman features largely 
unaffected. The mean effective temperature is T^ff - llOOOK 
for EF Eri, 12000 K for BL Hyi, and 10000 K for CP Tuc, with 
an estimated accuracy of about ~ 1500 K. Using this mean 
temperature in the tomographic analysis affects the theoreti- 
cal Zeeman absorption features only minimally, because their 
equivalent width reaches a maximum around 11000 K and 
varies little with effective temperature around the maximum. 
In passing we note that our mean effective temperatures con- 
firm the rather low temperatures of the white dwarfs in polars 
(Araujo-Betancor et al. 2005, and references therein). 

For the tomographic analysis, the flux and polarization 
spectra were collected into n = 4 phase bins for EF Eri and 
n - 5 for BL Hyi about equally spaced to cover the whole or- 
bit. The spectra from the half orbit of CP Tuc were gathered 
into three bins. 

3. General approach 

We determine the global surface magnetic field structure us- 
ing the Zeeman tomographic procedure described by Euchner 
et al. (2002, 2005, 2006). The process involves the inversion 
of the one-dimensional time series of rotational phase-resolved 
Zeeman flux and circular polarization spectra obtained in a 
non-accreting (low) state into a two-dimensional field distri- 
bution over the surface of the star Because of the finite signal- 
to-noise ratio, this inversion problem may allow more than one 
solution within the observational uncertainties (for a discussion 
see Euchner et al. 2002). This ambiguity arises from the fact 
that different models may have similar frequency distributions 
of the field strength B and the angle tj/ between field vector and 
line of sight and, hence, yield similar Zeeman spectra, but differ 
in the arrangement of the field over the surface. In this situation, 
it is advantageous to use the field vector in the accretion spot as 



deduced from cyclotron spectroscopy and broad-band circular 

polarimetry in a high state as a fixpoint and effective constraint 
on the tomographic procedure. In the present paper, we have 
not included such a constraint in a formal way, but use it to 
select between solutions of the tomographic process obtained 
for different assumptions on the field geometry. To this end, we 
follow the surface field outward and determine the maximal ra- 
dial distance reached by each field line. Field lines extending 
to more than - 107?wd are considered 'open'. The accre- 
tion stream can couple to field lines which reach out sufficiently 
far, with the actual radius at which coupling can occur depend- 
ing on the ram pressure of the accreting matter and the local 
field strength. The requirement that field Unes which originate 
at a specific point at the surface reach out to more than several 
white dwarf radii can effectively discriminate between different 
field models. To be sure, a model which provides a good fit to 
the Zeeman spectra and possesses field lines reaching far out 
at the required position may not be the correct model, but is 
as close to reality as we can presently get. As a further caveat 
note that the actual coupling conditions have not been investi- 
gated and part of the far-reaching field lines may not be acces- 
sible to the stream. Consequently, only a fraction of the long 
ribbon-like structures of far-reaching field lines which appear 
in some models may act as foot points of accreting field lines. 
Nevertheless, with this information included, the analysis of ac- 
creting white dwarfs may yield more definite results than that 
of isolated white dwarfs. 

As in our previous papers on isolated white dwarfs 
(Euchner et al. 2005, 2006), we fit the data with either a hy- 
brid model or a multipole expansion truncated at a maximum 
degree / = /max including all /max(/max + 2) components with 
m - ... Z. As a hybrid model, we consider the superposi- 
tion of zonal (m - 0) multipole components which are allowed 
to be inclined to each other and to be offset from the center 
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Table 2. Best-fit magnetic parameters for the truncated multipole expansions up to degree /max = 5 for EF Eri and up to /max = 3 for BL Hyi and 
CP Tuc. The coefficients and are in MG. The tilt angle of the multipole axis relative to the rotational axis is 74°, 32°, and 23° for EF Eri, 
BL Hyi, and CP Tuc, respectively. 
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of the white dwarf. Examples are, e.g., an offset dipole or the 
sum of dipole and quadrupole etc. Such combinations can eas- 
ily be visualized given the polar field strengths and orienta- 
tions of the individual components. In case of the multipole ex- 
pansion, on the other hand, the parameters of the basic dipole 
are easily interpreted, but the structure created by the higher 
multipole components is more difficult to judge (Euchner et al. 
2002). The hybrid models correspond to special situations not 
encountered in truncated multipole expansions of low /„,ax and 
we lack information from dynamo theory on the feasibility of 
such models. For the sake of economy and simphcity of pre- 
sentation, we present results for the offset dipole as a simple 
and popular model and for the fuU multipole expansion for ei- 
ther /max = 3 or 5, with occasional comments on other models 
(a multipole expansion up to /max = 4 was not tested). We use 
two graphic forms to present the results: (i) the 'fi-t^ diagrams' 
which depict the frequency distribution of field vectors over the 
surface of the star at a given orbital phase in the B-cos plane; 
and (ii) actual images of the field distribution. The latter in- 
clude (a) the field strength B, (b) cos t^^ = B;/B with 5/ the 
field component along the line of sight, and (c) an image of 
the maximum radial distance to which a field line extends that 
originates from a certain location on the star. 

We subjected the flux and polarization spectra at all orbital 
phases simultaneously to the tomographic analysis, weighing 
all wavelengths equally except for narrow intervals around the 
Balmer emission lines. An improved fit can be obtained by re- 
stricting it to the set of flux and polarization spectra at a single 
phase, but if the model parameters deduced for different phases 
disagree there is no unique solution (Euchner et al. 2006). 



4. Results 

The models are fitted to the average flux and circular polariza- 
tion spectra in four orbital phase bins for EF Eri, five bins for 
BL Hyi and three bins for the half orbit of CP Tuc. The orbital 
phase conventions used in this paper are the dip ephemeris for 
EF Eri (Piirola et al. 1987a) shghtly updated by including the 



ROSAT PSPC dip timings from July 1990 ' (Beuermann et al. 
1 99 1 ), the ephemeris for the start of the bright phase for BL Hyi 
(Wolff et al. 1999), and the dip ephemeris for CP Tuc (Ramsay 
etal. 1999). 

In our previous papers on the field structure of single white 
dwarfs, we considered the inclination of the line of sight rel- 
ative to the rotation axis as a free parameter of the fit. For 
the mCVs, however, independent and better information on / 
is available from the light curve and broad-band polarization 
studies in their high states. We use / = 55° for EF Eri (Piirola 
et al. 1987b; Achilleos et al. 1992) and / = 40° for CP Tuc 
(Thomas & Reinsch 1996; Ramsay et al. 1999). For BL Hyi, 
we use / = 32° (Schwope et al. 1995). 

Our Zeeman tomography uses the observed flux and circu- 
lar polarization spectra of all orbital phases. For conciseness, 
however, we show the spectra only for two selected phases and 
the globes representing the field structure only for ^ when 
the main accretion spot most directly faces the observer, i.e. 
00 - 0.10 for EFEri, (p^ ^ 0.20 for BLHyi, and probably 
00 - 0.50 for CP Tuc. Our phases closest to 0o are = 0.14 
for EF Eri, = 0.19 for BL Hyi, and = 0.43 for CP Tuc. 

' To = HJD244 3944.9518(6) -I- 0.056265949(14) E, with 90% con- 
fidence errors. 



Table 3. Dipole field strength fijjp, most frequent photospheric field 
strength fipmb, and ranges Smin to Bmax in the models of EF Eri, 
BL Hyi, and CP Tuc in MG. The former refers to the orbital phase 
when the accretion spot faces the observer (Fig. 3-5), the latter to all 
orbital phases combined. 



Object 


Model 


Bdip 


^prob 


^miii 


^max 


EFEri 


offset dipole 


44.0 


12 


10 


110 




'max = 5 multipole 




13 


4 


111 


BLHyi 


offset dipole 


59.5 


17 


15 


88 




/max = 3 multipole 




18 


13 


110 


CP Tuc 


offset dipole 


19.8 


10 


9 


41 




/max = 3 multipole 




10 


1 


69 
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Fig. 3. Visualization of the three magnetic field models discussed for EF Eri. The results shown are for orbital phase (p = 0.14 near the peaks 
of the X-ray and infrared cyclotron fluxes and shortly past the infrared absorption dip ((p = 0). The models are (i) the offset dipole (top), (ii) 
the full multipole expansion truncated at /,^iix = 3 (center), and (iii) the same for /^ax = 5 (bottom). The left panels show the B - if/ diagrams in 
which the frequency distribution of field vectors over the surface is represented by a grey scale. Also shown are the weight distributions (we) 
projected onto the B and the cos if/ axes. The right panels display globes with the distributions of the absolute value of the field strength B, of 
cos if/ = Bi /B with 6/ the field component along the line of sight, and of the maximum radius in units of /?wd to which the field lines extend. 
Lines reaching 10i?wD were not followed further out. The color bars above the globes indicate the range of the respective parameters over the 
visible hemisphere (see the online version of the paper for colored figures). Progressing orbital phase corresponds to an anti-clockwise rotation 
of the globes. 



Figures 1 and 2 show the flux and circular polarization 
spectra (grey curves) for the value of <p closest to <po and for 
another phase selected to point out differences in the Zeeman 
spectra. The spectral data in Figs. 1 and 2 are shown twice 
along with the best fit theoretical spectra for two models, the 
shifted dipole (model A) and the multipole expansion with 
Im^x = 3 or 5 (model B). The ordinate scales refer to the bot- 
tom spectra, the other ones being shifted upwards by arbitrary 
amounts. Since the fit to Ha cr+ always mimics that of the o"" 
feature we have omitted the former in the figures to avoid ex- 
cessive compression in wavelength. That feature is included in 
the fits, however 

We judged the fits by eye and by a formal global for the 
flux and polarization spectra at all phases combined and with 



all wavelengths weighted equally. As discussed by Euchner et 
al. (2005) the formal reduced are large, because the adjust- 
ment of the observed continua to the model continua is not per- 
fect and the standard deviations used in calculating account 
for the statistical noise in the data but not for the remaining 
systematic differences between model and data. We quote the 
global xl values which serve as a guide line, but also judge 
the merits and failures of individual models by eye and de- 
scribe them in words. Not surprisingly, significantly reduced 
xl are obtained by excluding the poorly fitting spectral regions 
from the fit. However, since these regions differ from object to 
object, we have refrained from including such restriction in a 
general way. We have assured ourselves, however, that the in- 
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Fig. 4. Visualization of the two magnetic field models discussed for BL Hyi, the offset dipole (top) and the = 3 multipole expansion 
(bottom), both shown for orbital phase (p = 0A9 with the main accretion region on the lower hemisphere in front. The B - i/f diagrams are on 
the left and the model field distributions on the right. Further details are as in Fig. 2. 



elusion of the poorly fitting regions does not affect the selection 
of the best-fitting model as the one with the lowest xl- 

For all three objects, the dipole, even after allowing for an 
off-center shift in the three spatial coordinates, does not provide 
a good fit, whereas the multipole expansions fare decidedly bet- 
ter. While it is possible that substantially more complicated hy- 
brid models might be successful, we are limited in the number 
of models that could be tested by the slow convergence proper- 
ties of our code (see Sect. 5). We quote some of the parameters 
of the dipole models in the text and list the coefficients and 
/i"' of the best-fit multipole expansions (Euchner et al. 2006) 
in Tab. 2. These coefficients are given in MG and, although 
such models are difficult to visualize, the numbers allow some 
insight into the field structure: the three / - 1 coefficients com- 
bine to define the dipole which is allowed to be inclined rela- 
tive to the multipole axis; the l,m = 2,0 coefficient describes 
the quadrupole aligned along the multipole axis and the follow- 
ing m - coefficients the octupole and higher multipole zonal 
components; the m + (tesseral) components are modulated in 
azimuth in addition to zenith angle. 

Figures 3 to 5 contain the representations of the field struc- 
ture at (/> ^ (/>() for both the offset dipole and the multipole 
expansions. In these figures, progressing phase corresponds to 
anti-clockwise rotation and a motion of a feature on the globes 
from left to right. The most probable field strength seen at the 
face when the prospective accretion spot faces the observer 
and the range of field strength over the whole star are listed 
in Tab. 3. We discuss the results in these figures and tables with 
each object below. 



4A. EFEri 

The main accretion region in EF Eri faces the observer near 
phase (p - 0.10 and is located ~ 30° from the rotational 
axis (Beuermann et al. 1987; Piirola et al. 1987b). The field 
strength in the accretion spot is low as indicated by the fea- 
tureless optical cyclotron continuum. Zeeman halo absorption 
in the ordinary ray of the cyclotron continuum (Oestreicher 
et al. 1990) and cyclotron humps in the infrared (Ferrario et 
al. 1996; Harrison et al. 2004; Howell et al. 2006) suggest a 
field in the range of 10 to 21 MG with different studies favoring 
different values. The positive circular polarization in the high 
state (Piirola et al. 1987b; Oestreicher et al. 1990) implies 
that the accreting field line in the main (X-ray emitting) spot 
points away from the observer (negative cos i^). There is evi- 
dence from spectropolarimetry (Oestreicher et al. 1990) and 
from broad-band polarimetry (Piirola et al. 1987b) that a sec- 
ond accretion region is located near the same meridian at a 
colatitude of about 115° (lower hemisphere) with a disputed 
polarity. While Piirola et al. (1987b) argue for the same polar- 
ity as the main spot, Oestreicher et al. (1990) favor opposite 
polarity. Both accretion regions are only about 80° apart, an ob- 
servation that led Meggitt & Wickramasinghe (1989) to suggest 
a quadrupole field. 

The best-fit offset dipole (model A in Fig. 1 and top panels 
in Fig. 3) has a polar field strength (of the unshifted dipole) of 
44.0 MG, is inclined to the rotation axis by 75°, and is shifted 
off center by Xd,3'd, andzd equaling 0.26,-0.01, and-0.177?wD, 
respectively, with /?wd the white dwarf radius. The z-offset 
is directed along the dipole axis, x and y denote perpendicu- 
lar directions (Euchner et al. 2002). A moderately good fit at 
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= 0.14 it accompanied by an utter failure at = 0.89 with a 
global ;^'y= 19.97. This is seen in the Ha and Hj3 cr^ features in 
the flux spectrum of Fig. 1 at 4800 and 6100A and in the polar- 
ization spectrum at all wavelengths. As an added disadvantage, 
the model lacks open field lines at the position of the accretion 
spot near the meridian in the upper hemisphere at ^ = 0.14 
(Fig. 3, top panel). As expected, setting the shifts perpendicu- 
lar to the axis equal to zero does not improve the situation. A 
more complex hybrid model of ahgned dipole, quadrupole and 
octupole, centered or ofF-centered, did not fare well either. Our 
fits quite reUably exclude this class of models. 

Of the multipole expansions up to l^^^ - 3 (with 15 fit pa- 
rameters) and /max - 5 (with 35 fit parameters), the latter yields 
a shghtly better fit with 13.99 vs. 14.88 for flie former. The 
'max - 5 case is shown in Fig. 1 as model B. There are some 
subtle differences resulting in a better fit of either model to one 
or the other spectral feature and, in summary, one would have 
to conclude that there is no good basis for including the added 
parameters of the /max = 5 model if judgment is based solely 
on the flux and polarization spectra. In line with this, we find 
that the B - t// diagrams in Fig. 3 (center panel for /^ax - 3; 
bottom panel for /^ax - 5) are similar in the predominance 
of field strengths of 10-15 MG, but differ in the extensions to 
higher field strengths. The decisive difference of the two mul- 
tipole expansions is depicted in the rightmost globes in Fig. 3. 
The /max - 3 model possesses only a small spot from which 
field lines reach far out. Such field lines, whether they close 
within the Roche lobe of the white dwarf or not are needed for 
the accretion stream to be guided towards the white dwarf. The 
field in this tiny spot is directed outward (cos tj/ > 0), however, 
leading to negative polarization of the cyclotron emission from 
this potential accretion spot, while the observed polarization is 
positive (cosi// < 0). The /max = 5 model, for comparison, dis- 
plays a large arc-like region of ingoing field lines which cov- 
ers also the expected location of the primary accretion region 
~ 30° from the rotation axis and facing the observer shortly 
past = (Beuermann et al. 1987; Piirola et al. 1987b). 
This ribbon of open ingoing field lines with cos tff < winds 
around a region of quadrupole-Uke low-lying magnetic arcs in 
the lower right quadrant at (p = 0.14^. The decisive point in 
favor of the /max = 5 truncated multipole expansion among the 
models studied is the correct polarity of the open field Unes 
with cos t/f <0 and a positive sign of the resulting circular po- 
larization of the high-state cyclotron emission from the main 
accretion spot. 

As seen from Tab. 2, the dipole component of the multi- 
pole expansion is relatively weak with a polar field strength of 
8.5 MG obtained by squaring the three dipole coefficients. The 
strongest components are azimuthally modulated quadrupole- 
Uke ones. All individual / = 3 - 5 coefficients are smaller 
than 10 MG, but their combined effect is significant in shap- 
ing the field structure (Fig. 3) which is not that of an m = 
quadrupole. Hence, from the present study, we conclude that 
the field structure of EF Eri is substantiaUy more complex than 



^ The interested reader may look up Fig. 7 of Schwope et al. (1995) 
which shows the field structure of such a region at its upper pole. 



that of a centered or offset dipole or quadrupole and may be 
even more complex than suggested by the present best fit. 

4.2. BL Hyi 

BL Hyi is another polar that displays a complex accretion ge- 
ometry, usually referred to as 'one-pole' and 'two-pole' ac- 
cretion in states of low and high accretion rates, respectively 
(Beuermann & Schwope 1989). The main hard X-ray and cy- 
clotron emitting accretion spot is located in the lower hemi- 
sphere of the white dwarf and is visible only for part of the or- 
bit. Its appearance at the limb of the white dwarf defines photo- 
metric phase = 0. It slowly disappears some 0.40 phase units 
later (Piirola et al. 1987a). The second emission region in the 
upper hemisphere facing the observer emits in the soft X-ray 
and EUV regime and is visible over much of the orbital pe- 
riod (Schwope & Beuermann 1993; Szkody et al. 1997). The 
sign of the circular polarization of the cyclotron emission from 
the main (second) pole is negative (positive) implying positive 
(negative) cos^/r (Cropper 1987; Bailey 1988; Schwope & 
Beuermann 1989). The dominant photospheric field strength 
is 22 MG, but the presence of significantly higher fields is in- 
ferred from Zeeman spectroscopy. A field of only 12MG was 
at times detected by Ha Zeeman absorption in the cyclotron 
continuum emission of the main accretion spot (Schwope et 
al. 1995), but the feature may have originated at some height 
above the white dwarf surface. 

Our present data cover the whole binary orbit in five about 
equally spaced intervals. The best-fit offset dipole model fea- 
tures a polar field strength of 59.5 MG, an inclination against 
the rotation axis of 80°, and offsets of -0.18, 0.15, and -0.03 
/?WD in X, y, and z. As in the case of EF Eri, the offset dipole 
model does not fare well. We find;jfj=200.3 which is so large 
because of the small standard deviation in the low-noise ob- 
served spectrum and obviously poor fits in a few places. The 
model fails to reproduce the Ha cr^ feature in the flux spectra 
near 6000A and also much of the detail in the circular polar- 
ization spectra over the entire wavelength range (Fig. 2, left 
panel). The fit to the polarization data is particularly poor at 
(p = 0.40 (not shown). We discard the offset dipole model be- 
cause it provides a comparatively poor fit and because of the 
lack of open field lines in both the upper and the lower hemi- 
sphere at(p = 0.19 at positions which might correspond to the 
main and the secondary accretion regions (Fig. 4, upper right). 

The multipole expansion truncated at /max = 3 fits better 
than the offset dipole, although there are still some systematic 
deviation in the flux and circular polarization spectra, in par- 
ticular near 4300A as well as 6050A and around 5500A, re- 
spectively. In spite of the still very high ;t'y=193.51, the field 
distribution looks promising. On the positive side, we note the 
correct representation of the reversal of the poliirization neiu" 
4600A between = 0.19 and <p = 0.81 and the approximately 
correct description of the Ha cr~ polarization near 6100A. As 
in case of EF Eri, the quadrupole-like components are strongest 
but, again, the other components are essential for the field struc- 
ture (Tab. 2). There is a longish ribbon of field lines reaching 
out moderately far (about 3/?wd)- It faces the observer near 
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Fig. 5. Visualization of the two magnetic field models discussed for CP Tuc, the offset dipole (top) and the = 3 multipole expansion 
(bottom), both shown for orbital phase (p = 0.43 with the main accretion region on the lower hemisphere in front. The B - i/f diagrams are on 
the left and the model field distributions on the right. Further details are as in Fig. 2. 



(p - 0.19 and yields negative circularly polarized cyclotron 
emission (positive cos i//) as required for the main accretion re- 
gion. A second region with negative cos i// may allow access as 
close as 35° from the rotational pole and may be responsible for 
the intense flaring soft X-ray emission (Schwope & Beuermann 
1989) and the positive circularly polarized cyclotron emission 
at orbital phases when the main spot is behind the white dwarf. 
The B - if/ diagram prominently shows field strengths between 
13 and 45 MG with a faint (and ill-defined) extension to beyond 
100 MG. The ribbon which may contain the main accretion 
spot crosses the 1 3 MG field minimum quite consistent with the 
featureless cyclotron continuum and the 12-MG Zeeman ab- 
sorption in the cyclotron emission of the main spot. These facts 
suggest that the model approaches reality, although the laig^xl 
lets us suspect that the final model is still different. The expe- 
rience from EF Eri with substantially different structures of the 
best fit multipole expansions for /max = 3 and /^ax = 5 suggests 
that such an improved model can be found. In this case, how- 
ever, convergence problems have prevented the construction of 
an Zmax = 4 or 5 model. We do not expect such a better-fitting 
model to possess a simpler structure than the /^ax = 3 model 
and conclude that BL Hyi possesses a complex field geometry, 
probably not completely unraveled with the model presented 
here. 

4.3. CP Tuc 

Unlike the majority of polars, CP Tuc was discovered by its 
hard X-ray emission (Misaki et al. 1996). It displays a broad 
dip in the X-ray flux which was used to derive a rotational 
ephemeris of the white dwarf (Ramsay et al. 1999). The narrow 



emission Une from the heated face of the secondary star yields 
the orbital period (Thomas & Reinsch 1996) which agrees with 
the rotational period proving synchronism. The ephemeris of 
these authors shows that inferior conjunction of the secondary 
star occurs at dip phase = 0.08. There is debate about the lo- 
cation of the accretion region: Misaki et al. (1996) suggested 
that the energy-dependent X-ray dip arises from photoabsorp- 
tion; while Ramsay et al. (1999) assign it to a self-eclipse 
by the white dwarf. In the first case, the accretion spot conve- 
niently faces the observer shortly before inferior conjunction; 
in the latter it faces away from the secondary. 

As noted above, our data cover only one half of the binary 
orbit. They were combined into two sets of flux and polariza- 
tion spectra at dip phases - 0.99 and 0.21 plus a noisier 
single set at <p - 0.43. In the Misaki et al. and Ramsay et al. 
interpretations, the accretion spot faces the observer near (p -0 
and (f) ^ 0.55, respectively. We show our spectra at = 0.99 
and (p = 0.43 in Fig. 2 (right panel). Model A is the offset 
dipole and model B the multipole expansion with /max = 3 of 
which the latter yields a formally better fit with;i'^ = 12.35 and 
11.78, respectively. Visual inspection, however, shows that the 
differences are not pronounced and that both models reproduce 
the major Zeeman features of H/3 and the higher Balmer lines 
in the flux and the polarization spectra reasonably well. The 
B-ip diagrams in Fig. 5 demonstrate the preponderance of field 
strengths close to lOMG which make CP Tuc a low-field polar 
and may explain the hard X-ray spectrum. The offset dipole 
has a polar field strength of 19.8 MG, is practically aligned 
with the rotation axis, and offset mainly in z by 0.21 7?wd- It 
has its high-field pole in the upper hemisphere, whereas the 
multipole model possesses low fields in the same region. The 
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dipole component of the multipole expansion has a polar field 
strength of 17.8 MG, but the quadrupole-like components are 
of similar strength and the octupole-Uke components are non- 
negligible (Tab. 2). Both models, offset dipole and multipole 
expansion, differ in the sign of the longitudinal field component 
over the near and far hemisphere. The featureless cyclotron 
continuum indicates a low field strength (Thomas & Reinsch 
1996) and the negative circular polarization of the continuum 
(Ramsay et al. 1999) shows that the accretion region must 
lie in the lower (upper) hemisphere for the dipole (multipole) 
model. Only the multipole model, however, possesses an ex- 
tended region of outgoing field Unes which faces the observer 
near = 0.55 as expected from the Ramsay et al. model (some 
0.12 in phase later than in the lower right globe of Fig. 5). In 
summary, CP Tuc is the third of the polars studied here with a 
field structure more complex than a simple ofi'set dipole. As a 
caveat we recall the incomplete phase coverage. The results for 
CP Tuc are, therefore, preliminary. 

4.3.1 . Summary of results 

For all three objects, the truncated multipole expansions yield 
significantly better fits than the offset dipole models. Only the 
former provide access to the surface of the white dwarf via 
field lines reaching sufficiently far out at the expected posi- 
tions on the surface. As seen from Tab. 3, dipole and multi- 
pole models have practically the same most frequent values of 
the field strength. They cover also similar ranges, in particu- 
lar, if one considers that the faint extensions in the probabiltiy 
distribution to the lowest and highest field strengths are not 
well defined. This similarity reflects the fact that both models 
more or less fit the principal features of the Zeeman spectra. 
Intuitively, one might consider rotating the offset dipole model 
to match the accretion conditions, but that changes the phase- 
dependent B - 4r diagrams and destroys the fits to the Zeeman 
spectra. Hence, the offset dipoles clearly cannot match all con- 
ditions simultaneously. In spite of the similar B-distributions, 
the field structures of the two models are significantly differ- 
ent and our conclusion in favor of stuctures more complex than 
offset dipoles is safe. 

5. Discussion 

We have presented the first Zeeman tomographic study of 
the field structure of white dwarfs in polars based on phase- 
resolved VLT spectropolarimetry. We have demonstrated that 
the studied stars possess field structures significantly more 
complex than simple centered or offset dipoles. Such a result 
was considered possible or even likely in many previous pub- 
lications, but detailed proof was not so far available. Our re- 
sults clearly demonstrate this complexity, although we have to 
caution that our best fits may not yet describe reality in every 
detail. The simplest parameter indicative of a structure more 
complex than a centered dipole is the range of field strengths 
over the surface of the star which exceeds seven for the best fit 
models for all three stars, while it would be two for a centered 
dipole. The presence of strong higher multipole components 
besides the dipole may be surprising considering the fact that 



the white dwarfs in CVs have a typical age in excess of 1 Gyr 
(Kolb & Baraffe 1999) sufficient to expect substantial decay of 
the higher order components (e.g. Gumming 2002). Recreation 
of higher order poloidal components from a toroidal interior 
field has been suggested by Muslimov et al. (1995), but it is 
not known whether the required strong toroidal field exists in 
magnetic white dwarfs. The single white dwarfs HE 1045-0908 
and PG 1015+014 (Euchner et al. 2005, 2006) have similarly 
complex magnetic field structures at an age of ~ 0.5 Gyr. At 
present, field evolution is not sufficiently constrained by obser- 
vations, but may become so when when the field structure of 
more objects becomes available. 

We find that fitting the field structure of accreting white 
dwarfs in CVs offers a decisive advantage over the analysis of 
isolated white dwarfs. The location of the accretion spot and the 
absolute value and direction of the field vector in the spot can 
be deduced independently, e.g., from X-ray and optical light 
curves and from broad-band polarimetry. Requiring that a suc- 
cessful field model compUes with this independent information 
turns out to be a powerful tool and is the main driver for our 
conclusion in favor of field structure more complex than an off- 
set dipole in all three objects. 

There is a semantic aspect worth mentioning. Comparison 
of the right-hand globes for the multipole expansions and the 
offset dipoles in Figs. 3-6 (including the phases omitted for 
conciseness) demonstrates that the latter have well-defined cir- 
cular regions of outgoing field lines which are properly ad- 
dressed as 'poles' . The corresponding regions in the multipole 
structures, on the other hand, are quite irregular and longish 
structures which imply accretion geometries which probably 
allow access to the white dwarf surface at more than one po- 
sition characterized by a wide range of angular separations. 
The resulting accretion geometry is no longer appropriately de- 
scribed by the dipole-inspired expressions 'one-pole accretion' 
or 'two-pole accretion' . 

Given the complicated field structures in the three polars 
studied here, it is desirable to extend the Zeeman tomographic 
analysis to a larger number of objects in order to distinguish 
between idiosyncrasies of individuals and the general proper- 
ties of the class. Such a program is feasible, but it calls also for 
a consideration of the limitations of our approach. One obvious 
limiting factor is telescope time. Although we have used typi- 
cally two orbital periods of high signal-to-noise spectra using 
the VLT and the spectropolarimetric capabilities of FORSl, the 
remaining noise in the circular polarization spectra limits the 
discrimination between the B - i/r diagrams of different models. 
A more extensive tomographic program requires to cover sev- 
eral orbital periods of each target with an 8-m class telescope. 
A second limiting factor is CPU time. In order to thoroughly 
test a given field model, typically 50 to 100 ;(f^-minimization 
runs are required becauce they tend to get stuck in secondary 
minima of the complicated landscape. Each step in this pro- 
cess requires to assemble the Zeeman spectra for the respective 
field model from the database. The resulting lack of speed is the 
main reason why we had to Umit the number of models tested. 

In the previous papers of this series (Euchner et al. 2005, 
2006), we have already discussed the alternative approach 
of Donati et al. (1994), the Zeeman Broadening Analysis 
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(ZEBRA). This method determines the best- fitting B - ^ di- 
agrams for each orbital phase interval directly from the data 
employing the Maximum Entropy Method (MEM) as a reg- 
ularization procedure. The advantage is predictably speed, the 
disadvantage is the uncertainty whether the individual B - i/r di- 
agrams are compatible with any global physical field model. 
The method faithfully reproduces the distribution in the abso- 
lute value of B over the visible hemisphere at a given orbital 
phase, but substantially smears the angle relative to the fine of 
sight (Donati et al. 1994). If interpreted in terms of a multi- 
pole model, the reconstructed B - i^-diagram then leads to spu- 
rious higher multipole components. Furthermore, the require- 
ment that the so derived diagrams contain the B - ^ combi- 
nation describing the accretion spot is easily implemented, but 
the exact location of the spot on the star is not, because the 
method forgoes imaging. On the other hand, testing a larger 
number of field models for compatibility with the B - tfr dia- 
grams for the individual phase intervals will probably be less 
time consuming than our present approach, because the large 
database of Zeeman spectra need no longer be accessed after 
these diagrams have been established. We plan to further study 
the difi'erent approaches in order to find the ultimately prefer- 
able one. 
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